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Abstract 

We present the phase diagram and dynamical correlation functions for the Holstein-Hubbard model at half filling 
and at zero temperature. The calculations are based on the Dynamical Mean Field Theory. The effective impurity 
model is solved using Exact Diagonalization and the Numerical Renormalization Group. Excluding long-range 
order, we find three different paramagnetic phases, metallic, bipolaronic and Mott insulating, depending on the 
Hubbard interaction U and the electron-phonon coupling g. We present the behaviour of the one-electron spectral 
functions and phonon spectra close to the metal insulator transitions. 
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Electron-phonon effects in strongly correlated elec- 
tron systems are expected to be significant but have so- 
far received little theoretical attention. We study these 
effects here using the Holstein-Hubbard model[l,2] 
where a local Einstein phonon mode couples linearily 
to local charge fluctuations. The Hamiltonian is given 
by 

H = Yl e ( fc ) c L c fc CT + U Y^ n ^ n n 

ka i 

+ cooY. b t b i +9Y. ^ + b * ) ( n *T + n 'i - X ) • 
i i 

We use a semi-elliptical band of width W = 4 and fo- 
cus on the particle-hole symmetric case at zero tem- 
perature. The phonon frequency is fixed to ujo = 0.2. 
We calculate the phase diagram, in the absence of long- 
range order, using a number of local approximations 
(DMFT-NRG/ED [2], DIA [3]). 

All methods agree on the phase diagram shown in 
the figure, ft consists of a metallic region surrounded 
by two distinct gapped phases, a Mott insulator for 
large U and a bipolaronic phase, when the electron- 
phonon coupling g dominates. The transition to the 



Mott insulating phase is always continuous, as indi- 
cated by the grey boundary curve; the critical value of 
U is largely independent of g. In contrast, the transi- 
tion to the bipolaronic phase is first order for U > 3, as 
indicated by the full black line, but is also continuous 
for smaller values of U. 

More of the physics of the interplay of the electron- 
phonon and electron-electron interactions is revealed 
in the spectra of the dynamic response functions. These 
were calculated using the DMFT-NRG [4]. Electron 
and phonon spectra are shown within the figure along 
different scans in the phase diagram. The results for the 
corresponding dynamic spin and charge susceptibilities 
are discussed in [1]. 

The two upper insets show the electron spec- 
tra pg(u}). The left-hand scan for U — 1 and g — 
0.0, 0.4, 0.46, 0.5 shows the continuous narrowing of 
the central resonance with increasing g and the sub- 
sequent opening of a gap at g c w 0.47. Contrasting 
behaviour is seen in the upper right inset (U — 5, g — 
0.0, 0.6, 0.7, 0.75) with an initial broadening of the 
central peak followed by its sudden disappearance at 
a critical g « 0.71 . 
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Fig. 1. Phase diagram for the particle-hole symmetric Holstcin-Hubbard model at zero temperature for a semiclliptical band of 
width W — 4 and a phonon frequency uo — 0.2. The insets show scans of the electron and phonon spectra, see text. 



The two insets below the electron spectra show 
the behaviour of the corresponding phonon spectra 
Pd(<*>) = — ■ilm((6 i ;6t^ . For U — 1 we see with in- 
creasing g a complete softening of the phonon peak as 
the transition to the bipolaronic phase is approached. 
The increase of negative spectral weight for ui < indi- 
cates a growing number excited phonons. At g c w 0.47 
we observe a two-peak structure whose low energy 
peak vanishes and whose high-energy peak hardens 
back to ujq after the continuous transition. For U — 5 
the softening of the phonons is much weaker due to 
suppression of charge fluctuations. There is no two- 
peak structure at the transition and again the peak 
hardens to ujq in the gapped phase. 

Similar features can also be seen in the scan for a 
fixed g — 0.45 and increasing U, shown in the right- 
hand inset. The transition from the bipolaronic to the 
metallic phase is reflected by a two-peak structure. 
Increasing U further suppresses charge fluctuations 
and hardens the completely softened phonon peak 
until the Mott transition is approached and the peak 



arrives at the bare frequency wo- 
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